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Abstract: A detailed theoretical investigation of polarization effects in superdeformed nuclei 
Q^ \ is performed. In the pure harmonic oscillator potential it is shown that when one particle (or 

hole) with the mass single-particle quadrupole moment q u is added to a superdeformed core, 
the change of the electric quadrupole moment can be parameterized as q e ff = e(bq u + a), 
and analytical expressions are derived for the two parameters, a and b. Simple numerical 
\ expressions for q e ff(q u ) are obtained in the more realistic modified oscillator model. It is 

also shown that quadrupole moments of nuclei with up to 10 particles removed from the su- 
perdeformed core of 152 Dy can be well described by simply subtracting effective quadrupole 
moments of the active single-particle states from the quadrupole moment of the core. Tools 
are given for estimating the quadrupole moment for possible configurations in the superde- 
. formed A ~ 150-region. 

i PACS: 21.60.-n, 21.10.Re, 21.10.Ky 
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1 Introduction 

o 

Superdeformed rotational bands were first identified |l[ from transition energies only. At an 
early stage, crude measurements of (electric) quadrupole moments were performed Q which 
were an important evidence for the large deformation. Recently, it has become possible 
to measure relative quadrupole moments with a much higher precision H, [|]. Since the 
total quadrupole moment depends sensitively on the specific orbitals that are occupied, these 
quadrupole moments have become an important tool to verify configurations of different 
bands. It is the aim of the present paper to study in detail how the quadrupole moment 
depends on occupied single-particle states in superdeformed nuclei. 

Adding a particle with a single-particle quadrupole moment, q u , to a core causes a change 
in deformation due to polarization effects. The size of the deformation change depends on 
the relative deformations, or quadrupole moments, of the core and the added particle. The 
deformation change can be translated to a change in the electric quadrupole moment and 
we shall refer to this change as the effective quadrupole moment q e ff- We will consider two 
different ways to calculate q e ff- One (microscopic) way is to calculate the total quadrupole 
moment from a sum of single-particle quadrupole moments at the minimum-energy deforma- 
tions for the configurations before and after the particle is added. The other (macroscopic) 
way is to calculate the quadrupole moment from a homogeneously charged body with the 
appropriate deformation and volume. The change of the electric quadrupole moment when 
one particle/hole is added then constitutes the effective electric quadrupole moment. 

Simple models Q suggest that a near-spherical Z = N nucleus changes its microscopic 
electric quadrupole moment with about q^/2 due to polarization, and by an additional q v if 
the added particle is a proton. Generalizing to a deformed nucleus we shall show that, in the 



pure harmonic oscillator (HO) model, q e ff can be written as 



(leff)p,n = e ( b P,n( e ) ■ Qu + a p ,n(e)) . (1) 

The parameters and a PiTl depend on deformation, e, and are different for protons (p) 
and neutrons (n). The relation (1) will be derived in section 2 and, for the microscopic 
as well as the macroscopic methods, analytic expressions will be given for b p ^ n and a P}Tl 
at superdeformation but also at other closed-shell configurations corresponding to different 
deformations. 

When comparing experimental and calculated quadrupole moments, the modified oscil- 
lator (MO) has often been used. In refs. ||, H|, the deformation for different superdeformed 
configurations was calculated using the Nilsson-Strutinsky cranking method and the quadru- 
pole moment was then obtained in the macroscopic way. The agreement between experiment 
and theory was generally found to be quite good. Furthermore, these quadrupole moments 
come close to those obtained in relativistic mean field calculations |7|]. 

In section 3 the quadrupole moment is studied using the MO potential. The polarization 
effects from adding one particle (or hole) to the yrast superdeformed 152 Dy configuration 
are studied in subsection 3.1 utilizing the macroscopic approach to calculate the quadrupole 
moment. The corresponding results from the microscopic approach are presented in sub- 
section 3.3. In both models we obtain simple relations corresponding to eq. (1), but with 
some modifications. The possible reasons for these modifications (cranking, hexadecapole 
deformation, Strutinsky renormalization, the P-term or the I ■ s-term in the potential) are 
analysed. 

In the models used here, superdeformed bands are understood from single-particle mo- 
tion in a rotating deformed potential. Then, it seems natural to ask if physical quantities 
can be described by adding effective contributions from particles in different orbitals. The 
relations found in ref. Q between high-iV^c configurations in superdeformed bands and the 
moments of inertia were based on the additivity of single-particle contributions. 
An attempt to test the additivity of experimental moment of inertia was performed in 
ref. |j. For the single-particle angular momentum, additivity of alignment was first tried 
in ref. [JOJ, then further tested in refs. [11, [l^] and found to work well. Similarly, it was 



concluded that specific orbitals lead to well-defined deformation changes corresponding to an 



additivity for deformations [11|. 



Deformations can be translated into quadrupole moments whose additivity were tested 
in selfconsistent Skyrme-Hartree-Fock calculations by Satula et al. Jk|. They extracted 
effective quadrupole moments from least square fits over a large number of configurations in 
the region of nuclei with Z = 64 — 67, N = 84—87, and found that the quadrupole moments of 
these configurations could be well described by summing the contributions from the orbitals 
involved. In the present paper, somewhat similar studies are described in subsections 3.2 and 
3.4 using the more simplistic cranked Nilsson-Strutinsky approach. However, our studies cover 
the whole A = 152 superdeformed region from 142 Sm to 152 Dy where a limited number of low- 
lying configurations are investigated in detail. Furthermore, comparisons are made between 
extrapolations based on calculated effective and bare single-particle quadrupole moments, 
respectively. Indeed, based on our analytical calculations in the HO, we get a microscopic 
understanding why additivity works so well also in realistic nuclear models. In subsection 3.5 
we compare with experimental data. Finally, a short summary is given in section 4. 



2 Quadrupole moments in the pure oscillator. 

The polarization effect on a deformed core by one particle is studied in the pure HO potential. 
There exist two equivalent methods to calculate this effect (see e.g. p. 510 in ref. pi): Either 



as a renormalization of the quadrupole operator due to the coupling of the single-particle 
excitations to the (isoscalar) giant quadrupole resonance (treated in RPA), or by considering 
the new, self-consistent deformation that results when one particle is added to a core. In the 
present paper, only the latter method is used. 

In subsection 2.1 some useful definitions are given, and an effort is made to calculate the 
polarization using double-stretched coordinates. In subsection 2.2 non-stretched coordinates 
are used and explicit expressions are derived for the polarization effect on an axially symmetric 
Z = N core. In subsection 2.3 the formulae are generalized to a core with Z ^ N. 



2.1 Double-stretched quadrupole moments for one kind of particles 

Consider a HO potential 



Vosc 



-m [uj x x + u> y y +u z z ) . (2) 



A specific orbital \v) = \n x n y n z ) is described by the number of quanta rij in the three 
Cartesian directions. The single-particle energies are given by e v and the single-particle mass 
quadrupole moments q u are calculated as 

q v = (v\ 2z 2 -x 2 -y 2 \ v). (3) 

In the potential, the equipotential surfaces are ellipsoidal with the axes proportional to l/o;?,. 
We will only consider axial symmetric solutions corresponding to cj x = u> y = u±. Elongation 
is then described by the standard parameter e fill , 

UJ Z = w II - y 

lo x = io y = iv (l + . (4) 
Volume conservation corresponds to 

u x u y uj z = (u} ) , (5) 

where the parameter loq is determined from the radius of the core. We now transform the 
physical coordinates (x, y, z) to dimensionless coordinates and furthermore introduce the 
stretched coordinate system , 



x = £ = W x, 

y m 

y' = v = \\-zrv, ( 6 ) 

V m 



z = C = A/ z, 

y m 

corresponding to the system where the eigensolutions separate in the three Cartesian direc- 
tions. The single-particle quadrupole moment becomes 

h 2 Z (2{n z + l/2) n x + l/2 n y + l/2\ ( _, 



moj \ l-2e/3 l + e/3 l + e/3 



For a system of Z particles (protons) in the potential, we calculate the total energy E 
as the sum of the single-particle energies under the constraint of volume conservation. The 
energy can be written as 

E = %oj x Yi x + hujyT,y + Twj z Yi z , (8) 
where the S's measure the total number of quanta in the different directions, 

s * = + 2) ; i = x >y> z - ( 9 ) 



The selfconsistent deformation, is obtained by minimizing the total energy, in the (e, 7) 
deformation space. The energy minimization is equivalent to a self-consistency between the 
potential and the matter distribution fj~4}| , namely that the ratio (x 2 ) : (y 2 ) : (z 2 ) is the same 
for the two distributions. This can be expressed as 

^x^x = ^y^y = ^z^z- (10) 

From now on we will assume axial symmetry corresponding to an equal number of quanta 
in the two perpendicular directions, 



Xy = rEi. (11) 



The axial symmetry corresponds to 7 = while e is obtained as (see e.g. refs. 17, |18| ) 

= 3(2E,-S ± ) 
6 4X z + £, ' 



(12) 



The total microscopic electric quadrupole moment Q mc is calculated as the sum of the 
single-particle quadrupole moments 

Q mic = e f q v = *1 A ( El ^ . (13) 

^Ti m w Vl-2e/3 l + e/3/ v ' 

The question which specifically interests us now is how Q mic is changed if one particle with 
a single-particle (mass) quadrupole moment q u is added to the core. Adding one particle to 
a spherically symmetrical HO potential (Yt x = T, y = Yi z \Q mtc = 0), it is well-known || that 
the total quadrupole moment becomes Q mic = 2eq v . 

Starting from an arbitrary axially symmetric deformation of the deformed HO, defined 
by the total number of quanta in the different directions, Ti z and Sj_, and the number of 
particles Z, a similar relation is found, however only when expressed in the double-stretched 
coordinates [19|, 

// hu± „ hu± „ hu z 

x = x, y = y, z = z. (14) 

m m m 

We define the single-particle quadrupole moment in these coordinates, 

qZ={v\2(z") 2 -(x") 2 -(y") 2 \u}. (15) 
It is then straightforward to show that at the self-consistent deformation eo defined by eq. ( |i~2|) , 

Q" = eY.l" = ^ (16) 
i.e. the matter distribution is always 'spherically symmetric' in the double-stretched system. 



When adding a particle with a double-stretched quadrupole moment q", we find a total 
double-stretched quadrupole moment, 

Q"(s ) = 2eq'l, (17) 

where Q" is calculated at the 'new' self-consistent deformation, but where £o indicates that the 
double-stretched coordinates are now defined with uj± and lo z corresponding to the original 
deformation, Eq, which is different from the new deformation obtained with the added particle. 
This is thus a generalization of the formula for spherical shape. 



However, it turns out that eq. (17) is not too useful when calculating how the physical 
quadrupole moment Q is influenced by the addition of a particle. It is straightforward to 
find linear relations between Q and (r 2 ) in the two systems (and also in the single-stretched 
system, eq. @), but the somewhat complicated form of these relations, and the fact that 
also (r 2 ) in the double-stretched system must then be analyzed, make us conclude that it 
is more straightforward to work directly with the formulae in the physical (non-stretched) 
coordinates. 

2.2 Quadrupole moments with a Z = N core. 

Consider a system of Z protons (with charge e) and N neutrons (with no charge) in their 
respective HO potential with the total number of quanta (see eqs. (pi 11)) described by 



) and (£ zn , E_i_ n ), respectively. The protons and neutrons are coupled in the stan- 
dard way used in the MO J2(|, namely that protons and neutrons have the same radius. For 
nuclear radii to be reproduced, this results in frequencies varying according to 
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D ( N-ZV' 3 D ( N-Z 



^•>=A^{^—) a A^{^— )• (18) 

where the standard value of D is 41 MeV corresponding to r Q = 1.2 fm. The rightmost 
expression has become standard in MO calculations and will be used by us except when 
discussing explicit Z and N dependences. We then note that 

1 ( N-Z\ x l' 3 _ (2Z) 1 / 3 



A 1 / 3 V A J A 2 / 3 

1 f l+ ^y< 3 _ (19) 



A 1 / 3 V A J A 2 / 3 

For the isolated systems of protons or neutrons, we can use eq. ( [12] ) to calculate the defor- 
mation of an arbitrary state, and thus the deformation change, 6e, if a particle is added. For 
simplicity we start by analysing a nucleus with an equal number of protons and neutrons in 
the core. By minimizing the total energy, i.e. the sum of the single-particle energies, of the 
Z = N proton-neutron system the equilibrium deformation is obtained as 

£ _l 3 (2T, zp — Sj_p) ^ 3 (2T, zn — £j_ n ) \ ^q-j 
2 \ 4S zp + 4S 2n + J 

This expression will be derived for the general Z ^ N HO system in subsect. 2.3 below, see 
eq. (pl|). Equation (|20|) gives the reasonable result that to lowest order in 5s, the addition of 
either a proton or a neutron leads to a deformation change Se/2, for the combined system. 
This is a simple model which can be studied analytically leading to closed formulae which 
should be helpful when considering shape changes and polarization effects in the more realistic 
MO model. 

Let us first note that the formulae in the previous subsection can easily be generalized 
to the Z = N system. Thus, if a particle with a mass quadrupole moment q v is added 



to a spherical system, because of its polarization is 'equally divided' between protons and 
neutrons, the total charge quadrupole moment will increase by 0.5eq v if the added particle is 
a neutron and by (1 + 0.5)eq u for a proton. Furthermore, eq. (|i~7| ) can be generalized in an 
analogous way for the Z = N system. 

For a deformed core we assume that the added particle has n z quanta in the axial direction 
and n x , n y quanta in the perpendicular directions. In an analogous way to the capital S's 
above (cf. eq. ©), describing the total number of quanta, we then introduce 

o z = n z + 1/2 

c± = n x + n y + l. (21) 
The single-particle quadrupole moment can now be written as, 

h 2 £ A 1 ' 3 ( 2a z a± \ 



(22) 



and the total microscopic electric quadrupole moment as, 



where the volume conservation factor is 




We can also define a macroscopic electric quadrupole moment Q mac calculated from an 
isotropic charge distribution with the total charge equal to Ze. For such a spheroid with 
the symmetry axis b and the perpendicular axis a, the quadrupole moment is given as 
(2/5)Ze(b 2 — a 2 ). Then with a radius parameter ro, i.e. a volume (4/3)7rrQ^4, and a quadru- 
pole deformation calculated from the S's according to eq. (p0|), the macroscopic electric 
quadrupole moment becomes, 



gmac _ LzeA 2/3 r% [ — ) | — ^ - | . (25) 

From eqs. (|25|,|23|), using eq. (Il9|) , we see that both Q mac and Q mtc are proportional to ZA 2 ^ 3 
as Xj_p and S 2p increase as Z^l 3 . Since e can be expressed in the S's, Q mac and Q mtc can 
be considered as functions of the independent variables Yi± and S z for protons and neutrons, 
respectively, together with the number of protons Z and neutrons N. 

We now determine how Q mtc (eq. (p3|)) and Q mac (eq. (|25|)) are changed from the addition 
of one proton by differentiating with respect to Z, T,± p and E zp , and from the addition of 
one neutron by differentiating with respect to N, Sj_ n and S 2n . The final expressions can be 
simplified quite a lot if instead of the original Ti± and S z (which are the same for neutrons 
and protons), we introduce the axis-ratio k between the z-axis and the perpendicular axes, 
which can be written as 

k=^~, (26) 

and use k and T*± as independent variables. For example, the volume conservation factor 
then takes the form 

^ = (27) 
uj 2k + 1 K ' 



We will use Aq for the number of particles in the reference nucleus which thus has Aq/2 
protons and Aq/2 neutrons. Then, with k measuring the deformation, q u takes the simple 
form 

qu = — A;" 1/3 ^o /3 (2fcn z + k + n z — l — N osc ), (28) 
Dm 

where instead of a z and a±, we have used n z and N osc = n x + n y + n z to characterize the 
particle. 

We will refer to the change in Q (Q mtc or Q mac ) when a particle is added as q e ff, which can 
generally be expressed in Aq, £j_ and k, in addition to n z and N osc , describing the properties 
of the added particle. These expressions are not very complicated but become even simpler 
if we put k = 2, corresponding to an axis ratio of 2 : 1, i.e. a superdeformed shape: 



iliac 



ff ) p = (A /2) 2/3 r 2 e(^1.6 + ^(1.2n 2 -0.3-0.6iV o 





(qTf7) n = (A /2y<* r z e [OA + ^ (1.2n z - 0.3 - 0.6iV osc ) 

( 'mic\ 



_L 

(A /2) 1/3 (8n z + 0.25 - 2.5iV osc ) (29) 

Dm 

-.2 



Note that these formulae are general for a superdeformed Z = N system in the sense that we 
do not require that the {Aq/2) lowest orbitals are filled; the only requirement is that £ 2 = T,± 
for the core. 

An interesting physical situation corresponds to the filling of the orbitals below the 2 : 1 
gaps of the HO. It is then important to note that only for every second gap, the deformation 
calculated (eq. (|l2|)) from the minimum of the sum of the single-particle energies corresponds 
to a 2 : 1 ratio of the nuclear axes (e = 0.6), i.e. it is only for these self consistent gaps that 
S 2 = Ej_. These are the Z = N = g = 4, 16, 40, 80, 140, . . . gaps, where we have introduced g 
for the particle number at these selfconsistent gaps. It is clear (see eqs. (29) and (p8[)) that 



if for one of these gaps, q e ff is plotted vs. q u with varying N osc for fixed n z , or with varying 
n z for fixed N osc , we obtain straight lines. If both n z and N osc are varied, the relation is not 
as evident, however. 

Of main physical interest are the orbitals close to the Fermi surface, and we note that 
the orbitals which are degenerate at 2 : 1 shape have a specific relation between n z and 
N osc , i.e. when N osc decreases by one, n z decreases by 2. Thus, for each of these bunches of 
degenerate orbitals, we will again get straight lines when q e ff is plotted vs. q v . The relation 
between N osc and n z is used to reduce the set of variables, and all the needed information 
about the number of quanta can be expressed by -/V^, which counts the number of shells at 
a deformation where the ratio uj z : uj± can be expressed by small integers. 

Before we write down the simplified expressions, we note that the relations at 2 : 1 defor- 
mation can easily be generalized to a k : 1 deformation where k is a small integer number. Self- 
consistent gaps are then formed for particle numbers Z = N = g = 2k, 8k, 20k, 40k, 70k, . . ., 
i.e. if the orbitals below these gaps are occupied, then S 2 = k^T<±. Furthermore, for degen- 
erate orbitals at k : 1 deformation, if N s h differs by (k — 1), then n z differs by k. For the 
selfconsistent gaps, g and Sj_ can be expressed in N^: 

9 = ^ (Nsh + 1) {N sh + 1 + k) (N sh + 1 + 2k) 



E± = i (N sh + 1) (N sh + l + kf (N sh + 1 + 2k) , (30) 



i.e. 

9 - 2k 

Si " (N sh + l + k)- (31) 

Using these relations together with eq. (|29| ) in its general form for an arbitrary fc-value, 
we obtain q e ff as functions in q u , eq. (|28|), for selfconsistent gaps at spherical (k = 1), 
superdeformed (k = 2), hyperdeformed (k = 3), etc. shape, 

2 xl/3 25 4/3 r 2 g r 2 e , 2g y/3 2 



^ ^-1) 0± 



/ mac\ _ 

V^/Jp 5S ± tflDrn 1 " 15\kJ v ~ 'V N sh + k + l 

(a™ c \ - pV /3 2g 4/3 rge r 2 e /2a\ 2 / 3 2 / 4 



fey 5£± h 2 /Drn 15\kJ v 'V + + 

n = L5e * -^(T) 1/3 ^(^ +t+i 4) < 32 > 

«1 - °- 5 ^ -^(t) 1/3 ^ (jv - + * +1±1) - 



All these functions are seen to be on the linear form 

q eff = e(bq u + a) (33) 

with 

a = a ± A (34) 

where ao is the average of q e ff for a particle and a hole with = 0, and A is the small 
deviation from this average due to particle or hole nature of the orbital, i.e. a particle added 
to one of the degenerate orbitals (with N osc = iV s /j + 1) just above the gap, or removed from 
one of the degenerate orbitals (with N osc = N s h) just below the gap. For prolate shape ao as 
well as A in eq. $34) are negative and the plus (minus) sign in this equation corresponds to 
a particle (hole). 

Equations ([32]) are valid for arbitrary k values. They are illustrated for the superdeformed 
(k = 2) g = 40 gap (N s h = 5) in Fig. 1. There we have primarily chosen to put the particle 
in an orbital above the gap with (N osc ,n z ) = (6,6), (5,4), etc. As (q^ffj has ao = (which 

holds for all selfconsistent gaps) A mac is directly seen as the deviation from for q v = 0. All 
a- values are negative, i.e. theg e j /-values are all negative for q u = 0. This is easily understood 
from the fact that the quadrupole moment of an added particle to a superdeformed core must 
be sufficiently large to induce an increased deformation. 

The slopes for the microscopic q e ff are the same as for the spherical case, i.e. (6 mic ) p = 1.5 
and (b mic ) n = 0.5. The slopes, i.e. the 6-values, in the macroscopic case depend on both k 
and g but, as illustrated in Fig. 2, in a way so that they are almost constant and close to one 
for all gaps (except for the very lowest ones). The 6 mac -values are the same e.g. for g = 40 
at 2 : 1 as for g = 20 at 1 : 1 (and g = 60 at 3 : 1). These are the third lowest selfconsistent 
gaps at each deformation. 

Using eq. ( |i~8|) with the standard value of D = 41 MeV an asymptotic value of 0.995 is 
found for b mac when N s /j — ► oo. It turns out, however, that the deviation from 1 is simply 
related to the value chosen for D. This value is generally determined so that R rrns = 1.2A 1 / 3 
fm , see e.g. [jTifl , leading to D = 41.2 MeV with one digit higher accuracy. Using this value 
of D instead gives b mac = 1.000 in the limit N s h — > oo. Furthermore, if when deducing 
D, we do not go to the asymptotic limit but require that R rm s = fm for each HO 

configuration, Tiujq will depend on g/k in such a way that b mac = 1 for all selfconsistent gaps 
at k : 1 deformation. Consequently, in these cases Fig. 2 just shows the effect of a constant 
Hloq for all particle numbers at each deformation. However, as D = 41 MeV is generally used 
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Figure 1: Effective electric quadrupole moments, q e ff, versus the single-particle mass quadru- 
pole moment, q u , at superdeformation, derived in the HO in the case when one particle is 
added to an orbital just above the gap. Solid and dashed lines are used for microscopic and 
macroscopic q e ff, respectively. Thick lines are used for the g = 40 gap and thin lines for the 
g = 80 gap (only drawn for neutrons). Very similar behaviour appears (not shown) when 
one particle is removed from an orbital just below the gap, i.e. A in eq. ( j34| ) is of minor 
importance. 



in MO calculations independent of particle number, we think it is interesting to illustrate 
the effect which this ^4-dependence leads to for the polarization of an added particle. In 
particular this means that, in microscopic calculations, R rms will be larger than 1.2A 1 / 3 fm 
for light nuclei. 

At spherical shape q e ff = for q v = 0, i.e. a = 0, as seen in eq. (|32|) (k = 1). For prolate 



shape (a m * c ) p n and (a mac ) n are negative. In Fig. 3 is shown how a varies with quadrupole 
deformation for a fixed number of particles. For 1 : 1, 2 : 1 and 3 : 1 deformation we have 
used the selfconsistent gaps while for 1 : 2, 2 : 3 and 3 : 2 deformation we have used the gaps 
whose deformation at minimal total energy is closest to the 'correct' deformation. For all 
deformations we have fitted a curve as a function of particle number (which becomes almost 
linear) to get the values at g = 40. 

In the high-iVsh limit and for selfconsistent gaps we get the asymptotic relations defining 
the macroscopic effective quadrupole moments, 

{b mac ) p = (b mac ) n = 0.995 

(a mac ) P = -(o±OA23^±g 1 A (35) 



(a mac ) n 



[ Q - 9U ^9 2/3 ± 0.423^^ f 



and the microscopic effective quadrupole moments, 

= 1.5 



{ h " dc ) r 
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Figure 2: Values of b mac (i.e. the slope of the dashed lines in Fig. 1.) as functions of number of 
particles of each kind (Z = N) in the nucleus. The stars connected by solid lines are the 1st, 
2nd, 3rd ... selfconsistent gaps for each of the deformations: spherical (1 : 1), superdeformed 
(2 : 1) and hyperdeformed (3 : 1). One can see that the values converge to an asymptotic 
value close to one (see text for details), and that it requires the same number of selfconsistent 
gaps (i.e. more particles at larger deformation) to get the asymptotic value. The circles 
show the values of the slopes for deformations (1 : 2, 2 : 3 and 3:2) where there are no 
selfconsistent gaps. Note that all values are quite close to the value at the asymptotic limit 
and that, as explained in the text, the differences from this value can largely be understood 
from the choice of the parameter hcoo. 

(°™*\ = -("■«3^/ / '±0.212^ 9 V3j (36) 
(«"")„ = -(o-306^ 2/3 ±0.212^ 9 V3). 



It is clearly seen that cto (eq. (|34|)) increases with g 2 / 3 (oc A 2 / 3 ) which is the same growth rate 
as for q u , and therefore the relative importance of ao is independent of particle-number. The 
polarization difference between particles and holes, A, on the other hand gets less important 
the heavier the nucleus is as it only increases as g 1 / 3 . 

Note that although the functions for macroscopic and microscopic q e ff are quite different, 
the sum of the q e f f values when adding one neutron and one proton is much more alike. This 
means that the total electric quadrupole moment for a nucleus with Z ~ N is approximately 
the same independently of if the microscopic or the macroscopic formula is used. 

The different behaviour of macroscopic and microscopic effective quadrupole moments 
have the same origin in the deformed case as the spherical case, and is caused by the two 
basically different assumptions applied. In the macroscopic approach, the protons as well as 
the neutrons are assumed to have a constant matter distribution inside the potential, i.e. also 
the matter distribution for protons and neutrons have identical deformations. Consequently, 
for a Z = N nucleus the addition of a proton or neutron has the same polarization effect, ex- 
cept for a constant factor caused by the added charge, i.e. same b- but different a-parameters, 
as seen in eq. ([^) for a deformation corresponding to a k : 1 prolate shape. Independent of 
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Figure 3: Dependence of a mac and a mic on deformation e. The values are given for g = 40 
but, as the values scale with the number of particles in the same way as the single-particle 
quadrupole moment, a/q u is (asymptotically) independent of the number of particles. The 
different behaviour for particles and holes (A in eq. (34)) has been ignored. 



deformation, we asymptotically obtain (b mac ) p = (b mac ) n « 1, see eq. ( |35[ ) (see also Fig. 1 for 
the k=2 case). 

To analyze the microscopic approach, we start from identical proton and neutron con- 
figurations corresponding to an equilibrium deformation e a . Then if a proton is added, the 
isolated proton system will get a new equilibrium deformation which we will refer to as e +5e. 
The deformation of the combined system is then e Q + 5e/2. The equilibrium deformation of 
the isolated neutron system is unchanged, e G . The electric quadrupole moment is then calcu- 
lated at the (non-selfconsistent) deformation e Q + Se/2. From eq, (12), it is easy to find out 
that in a pure proton system, half of the change in the quadrupole moment will be caused 
by the changed value of e, and half of the change by the change in the S's. In the present 
case with a proton added to a Z = N system, the change in the S's for the protons is the 
same as when a proton is added to a pure proton system while the change in deformation 
for the total system (the polarization) is half of that for the pure proton system. This gives 
(b mtc ) p = 1-5. If instead a neutron is added, there is no contribution from the change in the 
E's as they refer to the proton-configuration, but the deformation of the total system will 
still be e Q + 5e/2, and the contribution from this change in deformation will thus be the same. 
This gives (b mic ) n = 0.5, see eq. ( ^6|) (see also Fig. 1 for the k=2 case). 

The difference between the macroscopic and microscopic polarizations are thus caused 
by the lack of selfconsistency when the isolated proton and neutron systems have different 
equilibrium deformations. Then, in the macroscopic approach, it is assumed that the pro- 
ton and neutron matter distributions fully adopt the common deformation while they only 
do it 'halfway' in the microscopic approach. With the assumption that the proton-proton, 
neutron-neutron and proton-neutron attractions are the same, it seems that the microscopic 
approach should come close to a fully selfconsistent treatment. However, a stronger attrac- 
tion between unlike particles would correspond to a polarization in the direction suggested 
by the macroscopic approach. 

In subsections 3.1 and 3.3 we calculate q e ff in both methods for the modified oscillator 
potential, and in subsection 3.5 we compare the two methods with experimental data. First 



we will however study the polarization effects in the HO with a Z ^ N core. 



2.3 Quadrupole moments with a Z ^ N core. 

How will the results from the previous subsection change when the number of protons is far 
from the same as the number of neutrons? We will not derive all the formulae again but just 
investigate how different quantities depend on Z and N . To analyse quadrupole moments, 
the correct equilibrium deformation of the proton-neutron system has to be obtained. This is 
done by minimizing the total energy E(e, Z, N). The energy of the proton system or neutron 
system can be written as 



o 



1„, „ \ 1/3 



Ei = 3h uj 0i ( -H ±j H Z j ) ; j =p,n. (37) 



We know that %ojq follows eq. (|l9|), and £j_ and T, z are proportional to Z 4 / 3 and iV 4 / 3 for 
protons and neutrons, respectively. Altogether the energy becomes 

Z 5/3 

E p oc 



A 2/3 
AT 5 /3 



En oc (38) 

for the two different systems. In a region around their respective equilibrium deformation the 
proton and the neutron energies can be well approximated with parabolas. The total energy, 
which is just the sum of the two energies, can therefore be written 

Etot = E 0p + Cpie - Ep) 2 + E 0n + C n (e — s n ) 2 , (39) 

where the equilibrium deformation of each subsystem is 

^ _ 3(2S 2 j — 



4X^j + S^j 



j =p,n (40) 



and the energies at the minima (Eq p and Eo n ) and the stiffness parameters (C p and C n ), all 
with proportionality according to eq. (p8]), determine each parabola. By minimizing the total 
energy, Etot, the equilibrium deformation of the total system is obtained as 

Cp£p ~\~ C n £fi Z I Ep + N I £ n fA-\\ 

£ ° ~ Cp + C n ~ Z5/3 + AT5/3 • ^ 4i ^ 

The b values describe the part of q e ff which depends on single particle properties of the 
added particle (eq. (|33|)). They enter via the change 5e in equilibrium deformation, and 
for protons in the microscopic case by a direct contribution of 1 q v . The change in total 
quadrupole moment caused by the deformation change is to first order 

5Q def oc ZA 2 / 3 5e, (42) 

while the single-particle quadrupole moment is 

A 2/3 
A 2/3 

q vn oc j^f{a zn ,a± n ), (43) 



where f(a zp , <7_i_ p ) and f{a zn , crj_ n ) are independent of Z and N . In total we get the expressions 

y ^ x -7 42/3 Z 5 / 3 3{2a zp -a ±p ) 

I ft J p a ^/ 3 Z-V3 /(azp;fJ±p) - Z 5/3 +iV 5/3^^'^) 

where F(a zp , cr± p ) and F(a zn , a± n ) are independent of Z and iV. The previously derived 
6- values (eqs. ( ^B|j3^ )) can now be generalized to any proton- neutron system 

97-5/3 

(b™(Z,N)) p = z5/3 + N5/3 (b mac (Z = N )) p 
7 2/V 5 / 3 

(b mac (Z,N)) n = ^ z5/3 + jv5/3 (b mac (Z = N )) n 

{b miC (Z,N)) p = l+ z5/ 3^ 5/3 ((^o = iVo)) p -l) (45) 

keeping in mind that the Z = N 6-values were only deduced for selfconsistent gaps, but 
should be approximately valid for all Z- and N- values. For the nucleus 152 Dy the predictions 
from the HO is 6™ ac = 0.78, b™ ac = 0.93, 6™ c = 1.39, and b™ ic = 0.47. 

3 Quadrupole moments at superdeformation in the cranked 
MO potential. 

We will now continue to analyze polarization effects on quadrupole moments in the cranked 
MO model. Starting from the superdeformed 152 Dy configuration, we shall study how the 
quadrupole moment is affected by adding or removing protons or neutrons in specific orbitals. 
Contrary to the previous section, all calculations are performed at / ~ 40ft. It has previously 



been concluded [24] that at superdeformation, the general properties of the single-particle 
orbitals are almost unaffected by spin up to the highest observed values I = 60— 70ft. Further- 
more, our present study shows that the polarizing properties in the HO at superdeformation 
are essentially the same at / = Oft and at / = 40ft. In subsections 3.1 and 3.2 we study the 
macroscopic quadrupole moments, which have been used in recent comparisons with exper- 
iment and seem to work quite well j|, HJ, while similar calculations using the microscopic 
quadrupole moments are carried out in subsections 3.3 and 3.4. We perform a complete 



calculation in the cranked MO potential (parameters from ref. [21|) with Strutinsky renor- 
malization to get the proper equilibrium deformation, e and £4, of each configuration. The 
formalism of ref. |2^] is used which makes it straightforward to study specific configurations 
defined by the number of particles of signature a = 1/2 and a = —1/2, respectively, in the 
different N ro t shells of the rotating HO basis. The macroscopic electric quadrupole moment, 
Q mac , is obtained from an integration over the volume described by the nuclear potential 
and the microscopic electric quadrupole moment, Q mtc , by adding the contributions q v from 
the occupied proton orbitals at the proper deformation. One-particle polarization effects are 
investigated in subsections 3.1 and 3.3 while the additivity of several particles is studied in 
subsections 3.2 and 3.4. In subsection 3.5 the two theoretical methods are compared and 
confronted with experiment. 
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Figure 4: Changes in the total macroscopic quadrupole moment, q e ff, when a particle is 
added to or removed from the super deformed core plotted vs. the single-particle mass quadru- 
pole moment, q v . The circles and crosses are obtained from calculations in the Strutinsky- 
renormalized cranked MO with simultaneous energy-minimization in e and £4 direction using 
152 Dy as a core (values are given in Table 4). The quantum number N = N rot is indicated 
for each state. The solid lines, whose equations are given in the text, are linear least square 
fits to the respective points. The dashed lines are the quadratic least square fits for the HO 
potential (with £4 deformation included) with the {Z = 60, N = 80) core, as defined in Fig. 5. 



3.1 Macroscopic polarization effects of one-particle states 

In the present study, we first calculate the macroscopic quadrupole moments Q mac in neigh- 
bouring nuclei which only differ by one particle or one hole from the yrast superdeformed 
configuration in 152 Dy. From these Q mac -values effective electric quadrupole moments, q e ff, 
are obtained for different orbitals. For the orbitals around the Z = 66 and N = 86 gaps, 
q e jf is plotted vs. the mass single-particle quadrupole moment, q u , in Fig. 4 (see also Table 4 
below). It is evident that these values define approximate straight lines. By a least square 
fit we get the relation 

q eff = (1.02a, + 5.4) efm 2 (46) 

for protons and 

q eff = (1.30g„ - 43.7) efm 2 (47) 

for neutrons. This should be compared with 6™ ac = 0.78 and 6™ ac = 0.93 obtained from 
eqs. ( ^J45| ). It is evident that there are important differences. 

Let us analyze the differences between the Z = N HO and the MO results for 152 Dy by 
introducing, in the numerical calculations, the different terms one at the time. Starting with 
the Z = N HO we get the expected modifications when combining the HO superdeformed 
gaps Z = 60 and N = 80: Due to the neutron excess the average slope gets notably smaller 
than in the Z = N case and neutrons have a larger slope than protons, all in agreement with 
eq. (|45l). Cranking the system and performing a Strutinsky renormalization (i.e. including 
macroscopic liquid drop energy) only have a minor effect on the slopes, see dashed lines in 
Fig. 5. In the standard HO only quadrupole deformation e is used, while in the MO the 
energy is minimized in the quadrupole and hexadecapole deformation plane. In order to 
test the importance of the hexadecapole £4-deformation it is now included in the (cranked 
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Figure 5: Changes in the total macroscopic quadrupole moment, q e ff, when a particle is 
added to or removed from the Z = 60 and N = 80 core plotted vs. the single-particle mass 
quadrupole moment, q v . The curves are least square fits to the calculations in Strutinsky- 
renormalized cranked HO calculations. The dashed lines show the relations when the energy 
is minimized only in e deformation with £4 = 0. The solid lines (corresponding proton and 
neutron data marked with circles and crosses, respectively) show the relations when the 
energy is simultaneously minimized in the e and £4 directions. 



and Strutinsky-renormalized) HO, see the solid lines in Fig. 5. Due to the hexadecapole 
£4-deformation there is now a curvature in the relation q e ff = q e ff(q„), but the best linear 
fit is very close to the one obtained for quadrupole deformation only. 

When changing to the MO potential by including the / • s and ^ 2 -terms, there are large 
changes, see Fig. 4. The large differences between the HO and the MO are understood by 
considering the stiffness of the core in the two models. Around the equilibrium quadrupole 
deformation £0 of the core, its energy can be expressed as 

E core = E (e ) + C(e-e ) 2 , (48) 

and the energy of the added particle as 

e P art = e (eo) + K(e - e ). (49) 

The change in deformation of the total system, due to the added particle, is therefore K/2C. 
The stiffness C of the energy surface is larger for the HO than for the MO (see below), 
and this directly gives a larger deformation change in the MO for the same single particle 
quadrupole moment (same K). Consequently, in the MO the total quadrupole moment will 
increase more for high-iV osc orbitals with a positive deformation change and decrease more 
for \aw-N osc orbitals with negative deformation change. The slope of q e ff(qu) will therefore 
be considerably larger in the MO than in the HO. 

The net result of the change from the HO with a Z = N core to the cranked MO potential 
with the Z = 66, N = 86 core is thus an increase in the slope of the q e ff vs. q v relation for 
neutrons, but almost no change in the proton slope. Furthermore, the relations still become 
approximately linear also in the MO case. 

One might ask why the polarization effects are larger in the MO than in the HO. To inves- 
tigate this, we carried out calculations in the HO both with the HO closed shell configuration 



(Z = 60, iV = 80), and for the configuration corresponding to the yrast superdeformed band 
in 152 Dy. It turned out that, when corrected for the difference in mass, the different con- 
figurations had very similar stiffness. Calculations were then carried through for the 152 Dy 
configuration with only the Z 2 -term or the I ■ s-term included, indicating that the / 2 -term is 
responsible for approximately 60 % and the I ■ s-term for 40 % of the change in stiffness. 



3.2 Additivity of macroscopic effective quadrupole moments 

It is now interesting to investigate if the q e ff values are additive, i.e. if quadrupole moments 
of several-particle several-hole configurations relative to the superdeformed 152 Dy yrast state 
(with the Q mac value 1893 efm 2 ) can be calculated from the formula 

Qest = Q yDyyrast) + lef f ~ ^ Qef f ■ (50) 

particle hole 

The additivity will first be tested for excited states in 152 Dy where the number of excited 
particles are the same as the number of excited holes. We shall then test the additivity in 
other superdeformed nuclei all the way down to 142 Sm. 

Starting with 152 Dy, we give in Table 1 calculated deformations and quadrupole moments 
for a few n-particle n-hole configurations with rather low excitation energy. The quadrupole 
moments are calculated both by exact integration and by use of eq. (|50|) with q e ff taken 
from the 1-hole or 1-particle configurations given in Fig. 4. The agreement between these 
two methods is good with a typical difference of 2 efm 2 . 

Considering the good agreement between the 'exact' and 'estimated' macroscopic quadru- 
pole moments for 152 Dy, one could ask if similar methods could be used to estimate the 
quadrupole moment for other SD nuclei around 152 Dy. Furthermore, considering the approx- 
imately linear relation between q v and q e fj discussed above, it is interesting to investigate if 
these quadrupole moments can be obtained from a knowledge of Q for the reference nucleus 
and q u (but not q e ff) for the active orbitals. Quadrupole moments for configurations in the 
nuclei 150 Dy, 148 Dy, 150 Gd, 148 Gd, 144 Gd, 143 Eu and 142 Sm were therefore calculated in three 
different ways: 

• By direct calculation at the appropriate equilibrium deformation; Q eX act in Table 2. 

• From the quadrupole moment of the superdeformed 152 Dy yrast state and the sum of 
effective one-hole quadrupole moments; Q es t(q e ff) in Table 2. 

• From the quadrupole moment of the superdeformed 152 Dy yrast state and the sum of 
effective quadrupole moments calculated from single-particle quadrupole moments by 



the simple linear relations eqs. (46,47); Q es t{<lu) in Table 2. 



When choosing the configurations, we have to make sure that if two orbitals interact in 
the superdeformed region, both these orbitals should be either empty or occupied. This is the 
reason why the four N = 6 neutron orbitals (2 of each signature) are treated as one entity. 

In Table 2, Qest(leff) values are presented only for configurations which have holes in the 
orbitals used to get the effective one- hole quadrupole moments plotted in Fig. 4 while Q es t(Qu) 
are calculated for all configurations. The two estimates based on q e ff and q u , respectively, 
differ by up to 18 efm 2 for four particles removed. The difference between the calculated value, 
Q exact, and the single-particle estimate Q es t{<lu) is at most 34 efm 2 . The corresponding rms- 
value is 22 efm 2 for all configurations in Table 2. It is astonishing that the summation of 
effective quadrupole moments, calculated from eqs. (46, ^), describes the 'exactly' calculated 



values within 2% for 142 Sm, which is ten particles away from 152 Dy, and where the deformation 
has changed from e = 0.58 to e = 0.48 and e = 0.52 for the two studied configurations. 



Table 1: Calculated deformations and macroscopic quadrupole moments for SD configurations 
of 152 Dy. The two values of the quadrupole moment are obtained from a numerical integration 
(Qexact) and by adding and subtracting q e ff values to the yrast quadrupole moment (Q es t)- 
The values of q e ff are obtained from 1-particle or 1-hole configurations relative to the yrast 
state, see Fig. 4 and Table 4 below. 



Configurations of Dy. 


Qexact 

(efm 2 ) 


West 

(efm 2 ) 


e 


e 4 


yrast 


1893.2 




0.5820 


0.0166 


tt([651]3/2-)- 1 ([413]5/2-) 


1826.9 


1828 


0.5716 


0.0211 


7r([301]l/2 ) i ([532]3/2 ) 


1960.6 


1960 


0.6012 


0.0282 


z/([770]l/2+)- 1 ([402]5/2+) 


1749.1 


1750 


0.5519 


0.0135 


vr([651]3/2-)- 1 ([651]3/2+)- 1 
([413]5/2-)([413]5/2+) 


1765.2 


1766 


0.5608 


0.0241 


^([651]3/2-)- 1 ([651]3/2+)- 1 
([532]3/2~)([413]5/2-) 


1809.1 


1811 


0.5706 


0.0259 


tt([651]3/2-)- 1 ([413]5/2-) 
i/([770] 1/2+ ) - 1 ( [402]5/2- ) 


1686.4 


1685 


0.5415 


0.0182 


vr([651]3/2-)- 1 ([651]3/2+)- 1 
([532]3/2-)([413]5/2+) 
z/([770]l/2+)- 1 ([521]3/2-) 


1726.9 


1723 


0.5601 


0.0388 



The two 142 Sm configurations included have been measured in experiment |J and com- 
parison with these data will be discussed in subsection 3.5. In the Sm-bands, there are two 
orbitals which have to be handled with extra care. Thus, due to the difference in deformation 
the order of some orbitals closest to the Fermi-surface are different in 152 Dy and 142 Sm. At 
the deformation of 152 Dy, the proton N osc = 5 and neutron N osc = 4 closest to the Fermi- 
surface are 7r([532]5/2) and z/([404]9/2), i.e. the corresponding points in Fig. 4 are constructed 
from configurations with holes in these orbitals. At the deformation of the Sm-bands on the 
other hand, the 7r([541]l/2) and i/([411]l/2) orbitals are higher in energy and therefore, the 
bands are formed with holes in these orbitals relative to the 152 Dy bands. Consequently, their 
single-particle moments should be used in the relations to get the quadrupole moments in 
142 Sm, see Table 2. The orbitals are relatively pure in the 142 Sm configurations because, at 
the relevant rotational frequencies, the crossings between the orbitals occur at deformations 
somewhere between that of 152 Dy and 142 Sm. Therefore, these configurations are anyway a 
good test on how well additivity, based on q u , works 10 particles away from 152 Dy. 

3.3 Microscopic polarization effects of one-particle states 

In a similar way as in section 3.1 we now calculate the microscopic quadrupole moments 
Qmic _ e J2u=i If 111 neighbouring nuclei to superdeformed 152 Dy and plot q e jf vs. q u , see 
the solid lines in Fig. 6. Using the full MO, the relations no longer are approximatively linear 



Table 2: Macroscopic quadrupole moments calculated in three different ways (see text) for 
SD configurations in selected nuclei with one or several holes relative to the 152 Dy reference 
nucleus. 



nucleus 


configuration relative 


Q exact 


Qestiqeff) 


Qest(Qu) 


Q exact Qest(Qv) 




SD 152 Dy yrast 


(efm 2 ) 


(efm 2 ) 


(efm 2 ) 


(efm 2 ) 


152 Dy 




1893 








150 Dy 


v7~ 2 


1725 


1722 


1715 


10 


150 Gd 


tt6~ 2 


1735 


1733 


1722 


13 


148 Dy 




1701 


- 


1714 


-13 


148 Gd 




1576 


1562 


1544 


32 


144 Gd 


tt6- 2 u7~ 2 6- 4 


1381 




1365 


16 


143 Eu 


it6- 3 v7~ 2 6- 4 


1305 




1277 


28 


142 Sm 


vr6- 3 ([541]l/2)- 1 ^7- 2 6- 4 


1232 




1198 


34 


142 Sm 


tt6- 2 ([541]1/2)- 1 3- 1 
^7- 1 6- 4 ([411]l/2)- 1 










1419 




1410 


9 



but rather quadratic, 
for protons and 



q eff = (-0.008g 2 + 2.28g„ - 23.3)efm 2 (51) 



e// = (-0.009^ + 1.61g„ - 30.0)efm^ (52) 



for neutrons. This might suggest that it would be more proper to express q e ff not only as 
a function of the single-particle quadrupole moment, but also of the hexadecapole moment. 
However, as found below, the relations eqs. (^) and ( |52"D seem to work well in the limited 
region of superdeformed nuclei with A = 142 — 152, so at present we will make no attempt 
to generalize eqs. (51) and (|5z 



In a similar way as in the macroscopic case, the reason why the microscopic relations are 
so different from what was found in the Z = N HO calculations is now analyzed by performing 
calculations where the different terms in the full MO relative to the HO are introduced one 
at a time. 

First, allowing a different number of protons and neutrons and calculating also for non- 
selfconsistent gaps give linear relations with changes in the slopes in accordance with eq. (f45|). 
Introducing the Strutinsky renormalization and cranking only has a minor effect on the q e ff 
vs. q v relations. 

The result of including hexadecapole deformation £4 in the HO is shown in Fig. 7. As 
in the macroscopic case (see Fig. 5), the hexadecapole deformation introduces a curvature, 
which is, however, more than twice as large in this microscopic case. In Fig. 6 the MO result 
is compared to the HO result, both including the hexadecapole deformation. The two curves 
are seen to be rather similar, although the stronger polarization effect for the MO than for 
the HO, discussed above for the macroscopic case, can be seen also in the microscopic case. 
Furthermore, there is a small increase of the curvature of the q e jf vs. q u relation. 
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Figure 6: Changes in the total microscopic quadrupole moment, q e ff, when a particle is 
added to or removed from the super deformed core plotted vs. the single-particle mass quadru- 
pole moment, q v . The circles and crosses are obtained from calculations in the Strutinsky- 
renormalized cranked MO with simultaneous energy-minimization in the e and £4 directions 
using 152 Dy as a core (values are given in Table 4). The solid lines are quadratic least square 
fits to these points. The dashed lines are obtained from analogous fits for the HO potential 
with the [Z = 60, N = 80) core, as defined in Fig. 7. 

The reason why the introduction of the I • s and Z 2 -terms rather increase the curvature in 
the microscopic case (Fig. 6) but removes the curvature in the macroscopic case (Fig. 4) is 
not understood. 



3.4 Additivity of microscopic effective quadrupole moments 

The additivity in the microscopic case is checked, see Table 3, in the same way as in the 
macroscopic case, i.e. the electric quadrupole moment is calculated in three different ways 
Qexact, Qest(q e ff), and Q es t(qu), as described in subsection 3.2. 

The result of this test is that the additivity seems to work with a similar accuracy in this 
model as in the macroscopic approach. This is illustrated in Fig. 8 where the difference be- 
tween the calculated, Qexact, and the single particle estimated quadrupole moments, Q es t(qv), 
is plotted as a function of number of holes relative to the superdeformed 152 Dy core. 

For the three configurations where the quadrupole moment has been estimated based 
on both q e fj and q u , the maximum deviation is 7 efm 2 , see Table 3. The rms value for 
Qexact — Qest(qu) is 29 efm 2 for the eight configurations considered in Table 3, while the max- 
imum deviation is 53 efm 2 . In the microscopic case, summing effective quadrupole moments, 

), describes the two 142 Sm configurations with a 4 % accuracy. 



calculated from eqs. (|51 



3.5 



Comparison between macroscopic, microscopic and experimental ef- 
fective quadrupole moments 



and eqs. (51, ^)) between 



The two models studied give quite simple relations (eqs. (|4(| | 
q e jf and q v , but their predictions for specific orbitals are rather different as can be seen in 
Fig. 9. 
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Figure 7: Changes in the total microscopic quadrupole moment, q e ff, when a particle is 
added to or removed from the Z = 60 and N = 80 core plotted vs. the single-particle mass 
quadrupole moment, q u . The curves are least square fits to the calculations in Strutinsky- 
renormalized cranked HO calculations. The dashed lines show the relation when the energy 
is minimized only in e deformation with £4 = 0. The solid lines (corresponding proton and 
neutron data marked with circles and crosses, respectively) show the relation when the energy 
is simultaneously minimized in the e and £4 directions. 



In the macroscopic relations neutrons and protons with high q u - value have similar effects, 
while the microscopic relations give very different values for protons and neutrons, with 
exception for the lowest ^-values. In the microscopic case for neutrons the maximum q e ff 
is not obtained from the maximum q v . In Table 4 q u and Qeff~ values are given for several 
different orbitals close to the Fermi-surface in 152 Dy. The g e jj-values are shown both for 
the macroscopic and the microscopic case. For comparison values from Skyrme-Hartree-Fock 
calculations are also presented. By using these effective quadrupole moments together with 
eq. (|50|) (with the scalings discussed below) electric quadrupole moments can be estimated 
for a large number of superdeformed configurations in a rather large region around 152 Dy. 
There are also configurations with only q v - values given. From those one can estimate effective 
electric quadrupole moments by using eqs. ( ^J47| ) and eqs. (51,|52|), and through eq. (50) get 
good estimates of electric quadrupole moments for many more configurations in this region. 

In the comparison with experiments, see Table 5, the calculated quadrupole moments 
are scaled with factors to give approximately the same value for 152 Dy as the experimental 
data. This is partly motivated by the uncertainties in the absolute values obtained in ex- 
periment, due to the uncertainties in stopping powers. Also in Woods-Saxon, Hartree-Fock 
with Skyrme force, and cranked relativistic mean field calculations |2^, [l3|, 0] the values are 
systematically higher than in experiment. We see that the macroscopic method reproduces 
the relative changes with somewhat better accuracy in this region. The rms-values between 
the experimental and scaled theoretical quadrupole moments are 31 efm 2 and 48 efm 2 for 
macroscopic and microscopic models, respectively. The configuration with the largest dis- 
crepancy, v6~ 1 3~ 1 v7~ 1 in 149 Gd, has the largest discrepancy also in the Skyrme-Hartree-Fock 
calculation [13] and the error is almost the same. 

It is also interesting to note that the contribution to the change of quadrupole moments 
coming from protons and neutrons are very different in the two approaches. In the macro- 
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Figure 8: The difference between the calculated and estimated (based on q u ) quadrupole 
moment as a function of number of holes relative to the superdeformed 152 Dy core. The 
upper panel shows the macroscopic result from Table 2 while the lower panel shows the 
corresponding microscopic result from Table 3. 



scopic model 48 % of the total change in quadrupole moment from the superdeformed yrast 
band in 152 Dy to the exited band (7r6 _3 5 _1 ^7 _2 6 _ ) in 1 Sm comes from adding q e jf for the 
removed protons, while the corresponding number is 64 % in the microscopic model. 

In the mean field calculation based on e.g. the modified oscillator or Woods-Saxon poten- 
tial, it is required that the proton and neutron deformation are exactly the same, namely the 
proton and neutron single-particle potentials are defined for the same deformation param- 
eters. Then in the macroscopic method to calculate quadrupole moments, also the matter 
distribution is assumed to have the same deformation. As discussed in subsect. 2.2 for the 
HO this explains why for a Z = N nucleus (b mac ) p = (b mac ) n ~ 1, while the microscopic 
approach leads to different deformations leading to (b mic ) = 1.5 and {b mac ) n = 0.5. Then, as 
seen in Fig. 9, these expressions are modified by different factors but general features are still 
the same in Nilsson-Strutinsky-cranking MO calculations for Z < N nuclei. The numbers in 
Table 5 do rather support the macroscopic formula, e.g. when comparing 152 Dy and 149 Gd 
configurations with configurations in 151 Dy and 148 Gd which differ by one N = 7 neutron. 
This N = 7 neutron appears to have a large influence on the measured quadrupole moments 
in somewhat closer agreement with the macroscopic than the microscopic calculations. On 
the other hand, measured quadrupole moments in 131 ' 132 Ce [25] indicate a very small polar- 
ization for an N = 6 neutron in this region, even smaller than suggested by our microscopic 
calculations. We can conclude that more experimental data with high accuracy combined 
with comparison with selfconsistent calculations are required to disentangle the polarization 
properties of protons and neutrons, respectively. 



4 Summary 

We have investigated the polarization effects of a particle on a well-deformed core in the 
harmonic oscillator (HO) potential as well as in the modified oscillator (MO) potential. Two 
different ways to calculate the quadrupole moment (and thereby the polarization effect) were 



Table 3: Microscopic quadrupole moments calculated for SD configurations in selected nuclei 
with one or several holes relative to the 152 Dy reference nucleus. 
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considered. In the microscopic approach the electric single-particle quadrupole moments are 
summed at the appropriate equilibrium deformations, while in the macroscopic approach the 
quadrupole moment is calculated by considering the nuclear charge as uniformly distributed 
over its volume, again at the appropriate equilibrium deformation. Averaging over protons 
and neutrons, the two models were found to give similar results even though the individual 
proton and neutron contributions turned out to be rather different. 

In the pure HO model, it was found for a, Z = N system, that the change of the electric 
quadrupole moment when a particle (or hole) is added, q e ff, can be described by a simple 
linear relation in the single-particle mass quadrupole moment, q u : q e ff = e(bq u + a). Ana- 
lytical expressions were derived for the deformation and mass dependence of the parameters 
a and b. It turned out that in the microscopic model, b = 1.5 for protons and b = 0.5 for 
neutrons while in the macroscopic model, b showed some variation but was close to one for all 
deformations and particle numbers, for both protons and neutrons. These differences were 
explained from the way the proton and neutron matter distributions are assumed to adjust to 
each other when the equilibrium deformation of the individual systems are different. Allowing 
Z 7^ N ', neutron excess led to a decrease of the 6-values, especially for protons. 

In the macroscopic case, a was essentially equal to zero for protons. In the other three 
cases, it was negative for prolate shapes and positive for oblate shapes, and scales with mass 
A approximately in the same way as the single-particle quadrupole moment, i.e. proportional 
to A 2 / 3 . The fact that the parameter a is positive for oblate shapes and negative for prolate 
shapes is easily understood. The quadrupole moment of the added particle must overcome 
some value in order to increase the core deformation, and this value is obviously positive for 
prolate shapes, negative for oblate shapes and zero for spherical shapes. 

In the Nilsson-Strutinsky cranking calculations we used the MO potential as the micro- 
scopic potential, and calculated effective quadrupole moments around the superdeformed core 
of 152 Dy. Both the macroscopic approach and the microscopic approach were used. From a 
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Figure 9: Effective electric quadrupole moments versus single-particle mass quadrupole mo- 
ments calculated in the MO for one particle/hole outside the superdeformed 152 Dy core. 
Dashed lines are used for the macroscopic and solid lines for the microscopic method. It is 
clear that for specific combinations of protons and neutrons removed from the 152 Dy core, 
the two methods could lead to rather different results. 



basic point of view the microscopic way of calculating quadrupole moments appears most rea- 
sonable. On the other hand the macroscopic approach has been used frequently in previous 
realistic calculations and been found to work well. In the macroscopic model numerical calcu- 
lation indicated that the linear relation between q e ff and q v were valid and the polarization 
factors were then numerically obtained as b p = 1.02, b n = 1.30, a p = 5.4 fm 2 and a n = —43.7 
fm 2 . These numbers were different from the factors deduced from the pure HO in the same 
macroscopic way for an Z = 66, N = 86 superdeformed nucleus (b p = 0.78 and b n = 0.93). 
The reason for this deviation was explained as being due to the decreased stiffness of the 
potential energy surface around the minimum for the MO. In the microscopic model there 
appeared to be a stronger dependence on hexadecapole deformation which led to a need for 
quadratic relations. 

Additivity of effective quadrupole moments in superdeformed nuclei was investigated and 
found to work surprisingly well. Adding q e ff- values, calculated from one-hole and one-particle 
states outside a superdeformed core of 152 Dy, quadrupole moments could be well described 
in an extensive region of superdeformed nuclei. Similar results using the Skyrme-Hartree- 
Fock method were previously obtained by Satula et al. [|l3[. Furthermore, using the simple 
relations for q e ff as a function of q u , quadrupole moments could be estimated in an even 
larger region using only the total quadrupole moment of the core 152 Dy together with q u - 
values for the active single-particle orbitals as input. For example, in the macroscopic case, 
the 10-hole configurations describing two observed superdeformed band in 142 Sm were both 
estimated within a 2 % accuracy relative to the values obtained from a full calculation for 
these bands. In the microscopic case the additivity worked with a somewhat smaller accuracy 
and we obtained a 4 % accuracy for the two 142 Sm bands. 

From the (bare) single-particle quadrupole moments given in Table 4 it should be possible 
to estimate total electric quadrupole moments with a reasonably accuracy for configurations 
in a quite extended region of superdeformed A ~ 150 nuclei. 

The experimental data are reproduced in a good way by the theoretical calculations with 



somewhat smaller discrepancies using the macroscopic method, see also refs. ||, [|, [|. 

The surprising accuracy of the additivity suggests the possibility of a shell-model type 
description of superdeformed nuclei, utilizing a superdeformed core and a valence space con- 
sisting of superdeformed one-particle (one-hole) states. 

We are grateful to A.V. Afanasjev for useful comments on this manuscript. I.R. and S.A. 
thank the Swedish National Research Council (NFR) for financial support. 
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Effective quadrupole moments calculated for orbitals in the SD A = 150 region. The macro- 
scopic and microscopic values, calculated in 152 Dy ± 1 particle, are the ones used when 
producing Fig. 4 and Fig. 6. For other orbitals only g^-values is given. They can be used 
together with the relations eqs. (46, 47) and eqs. (51, 52) to estimate macroscopic and mi- 
croscopic effective quadrupole moments. The Skyrme-Hartree-Fock calculations are from 
ref. [13]. 
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(a) The orbitals z/([642]5/2) and i/([651]l/2) are mixed in the Dy-region and are not pure in 
152 Dy. If none or both orbitals of the same parity is unoccupied the quadrupole moment can 
be correct calculated from these values, else extra care should be taken. The labels are valid 
for the Skyrme-Hartree-Fock values. 



Table 5: Macroscopic and microscopic quadrupole moments are compared with experiments. 
The holes are specified relative to the reference, superdeformed 152 Dy yrast configuration. 
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The experimental data are from a ) ref. ||, b ) ref. Q, c ) ref. ||. In the configurations marked 
' the hole is forced to the i/([411]l/2) -1 orbit in order not to mix with zv([404]9/2) _1 orbit 
while in the configuration marked e ) the hole naturally comes in the i/([411]l/2) _1 orbit. 



